ACOUSTIC FLUCTUATIONS DUE TO 
SHALLOW WATER THERMAL MICROSTRUCTURE 



James Burgess Hagen 



cv KNOX UBRW 

naVM. POSTGRADUATE St 

SS5^«. 



;hool 

93940 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESS 




ACOUSTIC FLUCTUATIONS DUE TO 
SHALLOW WATER THERMAL MICROSTRUCTURE 

by 

James Burgess Hagen 
September 197^ 



Thesis Advisor: 



E.B. Thornton 



Approved for public release; distribution unlimited. 



lfl6488 



b 




UNCLASSIFIED 

SECURITY CLASSIFICATION OF this PAGE Data Entarad) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4 . TITLE (and Submit) 

Acoustic Fluctuations Due tb Shallow 
Water Thermal Microstructure 


5. TYPE OF REPORT ft PERIOD COVEREO 

Master’s Thesis; 
September 1974 


«. PERFORMING ORG. REPORT NUMBER 


7. AUTHORf#; 

James Burgess Hagen 


•. CONTRACT OR GRANT NUM B ER(*) 


9. PERFORMING ORGANIZATION NAME ANO ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


10. program element, project, task 
AREA ft WORK UNIT NUMBERS 


It. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

September 1974 


13. number of pages 

154 


U. MONITORING AGENCY NAME ft ADDRESS (// dlllarant /root Controlling Oiilci) 

Naval Postgraduate School 
Monterey, California 9 39 h 0 


15. SECURITY CLASS, (ol Ihla tiport) 

Unclassified 


15 a. OECL ASSIFICATION/DOWNGRADING 
SCHEDULE 


16. DISTRIBUTION STATEMENT (o 1 thla Raport) 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (ol tha abatrmct ant at ad In Block 30, U dlllarant tram Raport) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Contlnua on ravaraa aida if nacaaamry mnd Idantlfy by block n\xmbme% 

Acoustic Fluctuations 

Thermal Microstructure 

Shallow Water Thermal Microstructure 


20. ABSTRACT (Contlnua on ravaraa aida It nacaaamry and Identify by block numbar) 

Simultaneous measurements of waves, orthogonal water 
particle velocities, temperature, conductiv/ity and acoustic 
phase and amplitude measurements at 20 and 40 kHz over 
horizontal and vertical, 2m paths were made? at the NUC 
oceanographic tower, San Diego in April 19774. 

Four stochastic models for acoustic propagation including 
geometric. Born, Rytov and Debye approximation models were 



DD 1 jan M 73 1473 EDITION OF 1 NOV 65 IS OBSOLETE 
(Page 1) S/N 0102-0U- 6601 | 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS PAGE (Whin Data Entarao) 




UNCLASSIFIED 



CkCUKlTY CLASSIFICATION OF THIS P*GEr*^#n Ent*r*d? 



(20. ABSTRACT continued) 

tested. The Krosil'nikov - Obhukov formulation of the Born 
approximation yielded reasonable results. Diffulty in 
establishing scales obviated good results from the Rytov and 
Debye models. 

Taylor thermal microscales from spatial correlation 
analyses determined vertically across the thermocline varied 
from 6.8 to 178 cm. The measured variance of the index of 
refraction ranged from ^.5x10“^ to 2.7xl0“3. 

Temporal stochastic analyses showed horizontal acoustic 
path losses were ten times vertical losses because of 
diffraction effects. The primary mechanisms for acoustic 
variability were internal and surface waves pumping the 
thermocline up and down. 



DD Form 1473 (BACK) 
, 1 Jan 73 

S/N 0102-014-G601 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS P AGEfPTt^n Dmt* Enfrmd) 



Acoustic Fluctuations Due to 
Shallow Water Thermal MicrostTuctTire 



by 



James Burgess ^agen 
Lieutenant, United States Navy 
B.S.E.E., University of Colorado, 1966 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OCEANOGRAPHY 



from the 

NAVAL POSTGRADUATE SCHOOL 
December 1974 



//O^ 



ABSTRACT 



D’./Clp Y KNCX II-'- r 
NAVAL Pc. 3 1 



WO.41 Ej 1 w y, 



CnL,, 






Simultaneous measurements of waves, orthogonal water 
particle velocities, temperature, conductivity and acoustic 
phase and amplitude measurements at 20 and 40 kHz over 
horizontal and vertical, 2m paths were made at the NUC 
oceanographic tower, San Diego in April 1974. 

Four stochastic models for acoustic propagation including 
geometric. Born, Rytov and Debye approximation models were 
tested. The Krosil'nikov - Obhukov formulation of the Born 
approximation yielded reasonable results. Difficulty in 
establishing scales obviated good results from the Rytov 
and Debye models. 

Taylor thermal microscales from spatial correlation 
analyses determined vertically across the thermocline varied 
from 6.8 to 178 cm. The measured variance of the index of 

, _ii 

refraction ranged from 4.5x10 to 2 . 7 x 10 . 

Temporal stochastic analyses showed horizontal acoustic 
path losses were ten times vertical losses because of 
diffraction effects. The primary mechanisms for acoustic 
variability were internal and surface waves pumping the 
thermocline up and down. 
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I. INTRODUCTION 



A. HISTORY 

A continuing program of in situ ocean microstructure 
research has been conducted at the United States Naval 
Postgraduate School (NPS) since 1971. Under this program, 
a series of three annual experiments dedicated to the study 
of the complex temporal and spatial relationships between 
multiple parameters measured in shallow water have been 
conducted by the Department of Oceanography. Each of these 
experiments has been restricted in scope to measurements 
inside a 9 m^ volume of water and to oceanographic phenomena 
which occur in the frequency range 0.001 to 2.5 Hz. 

B. OBJECTIVES 

In this thesis, initial analysis of certain parameters 
measured in the latest of these experiments is reported. 
Primary emphasis is given to stochastic temperature relation- 
ships affecting high frequency (20 kHz) CW sound transmission 
over a short path. The particular objectives of this 
experiment were to investigate the dependence of acoustic 
phase and amplitude fluctuations on the magnitude and 
characteristics of the ocean thermal structure and the 
intensity of the related turbulence. 
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II. THEORY 



The description of oceanic turbulence by stochastic 
techniques has followed developments in the atmosphere. 

One of the principal difficulties in application of the 
theoretical methods provided from atmospheric research has 
been determination of the range of applicability of certain 
turbulence models in terms of ocean parameters. Another 
has been the application of these models to sound propagation 
in a manner similar to the relationships developed for 
electromagnetic wave propagation in the atmosphere as 
summarized by Strohbehn [1968]. Historically, attempts to 
resolve these difficulties have been interconnected since 
their primary relevance is to acoustic prediction which 
motivated the research. 

A. STOCHASTIC MODELS FOR THE MEDIUM 

Three expected value techniques which have been applied 
to turbulence calculations in the ocean are 1) the spatial 
correlation (autocorrelation) function, 2) the structure 
function and 3) a "piecewise correlation" function. 
Application of each technique requires a statistically 
homogeneous medium but the order above indicates their 
relative applicability to increasingly less stringent 
stationarity requirements. The function form for each is 
listed below. 
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The autocovariance function, R (p): 

X 



R X (P) 


= E[T(x)T(x+p)] . Cl) 


The structure 


function, D (p) 


D X (P) 


= E{[T(x) - T(x+p)] 2 } . (2) 


The piecewise 


correlation function, R'(p) 

X 


R ; ( p> 


_ cn E{[T(x) - T(x+P')] 2 } 

" hL1 2"R (0) 

X 

D (p») 

■ E[1 - rrw ] • (3) 

X 



E indicates expected value and p is the spatial lag. 

For completely stationary processes, the structure function 
and the "piecewise correlation" function reduce to 



D X (P) 


= 2[R x (0) - R x (p)] , (*0 


and R’(p) 


= E[1 - D x (p)'/2R x (0)] = R X (P')/R X (0). 

(5) 



The stationarity difficulties often encountered in the 
use of R (p) are partly overcome by the use of D (p) because 

X X 

the large scale effects in the ocean are differentially 
filtered out [Black 19653 . The piecewise correlation function 
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can be applied to measurements which have only "piecewise 
stationarity . " This makes it the most powerful tool of the 
three techniques. "Piecewise stationarity" is defined as 
stationarity (in the ordinary sense) over a sufficient 
record length for which the structure function and variance 
can be calculated. This method was proposed by Whitmarsh, 
Skudrzyk and Urick [1957]. 

Various analytical equations have been used to approxi- 
mate the form of the auto covariance function used to describe 
the environment. Liebermann [1951] used a simple exponential 
form in his early experiments. Others have used a Gaussian 
form with some success [Stone and Mintzer 1953]. Use of the 
exponential form yields the best fit to experimental data 
above certain spatial lags, but it underpredicts the number 
of "small" turbulent blobs. This undesirable effect is 
eliminated by the use of the Gaussian form. 

Kolmogorov [19*11] predicted a 2/3 power law dependence 
for the structure function for turbulent flow in which there 
was a decoupling of the energy producing and energy dissipa- 
tion regions. Batchelor [1953] transformed this result and 
demonstrated that in wave number space energy spectra have 
a wave number to the -5/3 power law dependence. This was 
first verified convincingly by Grant et al. [1962] for ocean 
measurements and Pond et al. [1963] for atmospheric measure- 
ments. Ocean measurements do not always shov; the dependence 
since the ocean tends to be stratified much more intricately 
structured than the atmosphere which makes the recording of 
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continuous stationary data unlikely. However the -5/3 
power law depends not so much upon isotropy as upon the 
clear separation of scales of ' energy production and energy 
dissipation in which a cascade of energy from larger to 
smaller scales occurs. Whether or not such a form is 
ubiquitous in the ocean is not as significant as the presence 
of horizontal sheets which could cause multiple reflection 
layers. Very sharp, multiple thermal layers have been 
reported in the Arctic [Denner 1971], for example. 

Recently it has been shown that surface wave effects can 
be removed from ocean turbulence measurements with the result 
that purely turbulent flow remains and a -5/3 power relation- 
ship may be observed [Thornton and Boston 197*0. These 
findings were based upon measurements both in and around the 
thermocline and showed that anisotropy was not a critical 
difficulty, provided additional energy sinks or sources did 
not result (e.g. in the thermocline). 

Unfortunately, other energy contributions to the 
turbulence spectra do exist. Microstructure related 
turbulence from salt-fingering is one form [Gregg 1973 , 

Denner 1971]. Billow turbulence resulting from shear 
instabilities related to internal waves has also been 
reported by a number of investigators [Woods and Wiley 1972, 
Gregg 1973, Powell 197*0. Further, the existence of the 
internal wave itself could be expected to add energy at low 
frequencies as could surface wave interaction with the 



thermocline . 



B. WAVE PROPAGATION MODELS 

Four methods derived from the Helmholtz equation have 
been used to describe wave propagation through the turbulent 
ocean. These are: 1) geometrical optics, 2) the Born 

approximation [Born 1933], 3) Rytov's method [Rytov 1937], 

and 4) the Bebye approximation [Frisch 1968]. While it is 
not a necessary assumption, most applications of these 
methods presume that the ocean is homogeneous so that the 
scalar wave equation can be used. Such an assumption places 
a rather severe restriction on the acoustic wavelength, X, 
such that 



X << Z (6) 

where Z = inner scale of the turbulence. This restriction 
o 

will be discussed in detail later. 

1 . Geometrical Optics 

The geometrical optics (ray theory) approach is a 
linearization of the wave equation which assumes only 
Rayleigh scattering and refraction effects are important. 
Since diffraction effects are ignored, the high frequency 
(short wavelength) restriction, 

X << £ q 2 /L (7) 

where L = acoustic path length, is placed on the results. 
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One of the first experiments which used a stochastic 
technique to examine the acoustical effects of small-scale 
temperature inhomogeneities in the ocean was conducted by 
Liebermann [1951]. Through the use of the classical optical 
model, Liebermann was able to predict certain scintillations 
in sound amplitude observed at sea by the equation: 



V 


= /if a 2 (L/a)^ 

15 


(8) 


V I 2 


= E[(^) 2 ] 


(9) 


2 

a 


= E[(|^) 2 ] 


(10) 



o 



2 

where Vj = variance of the sound intensity, I = mean sound 

intensity, AI = the instantaneous deviation of sound intensity 

2 

from the mean, a = variance of the index of refraction, 

L = acoustic path length and a = Taylor temperature micro- 
scale [Lumley and Panofsky 1964], Liebermann also showed 
that sound reverberation encountered in the ocean was larger 
than that predicted due to thermal inhomogeneities alone. 

This was substantiated by later developments which have 
shown that biological scattering, chemical relaxation and 
bubbles in the near-surface region are important contributors 
to sound attenuation. 

Sagar [1955,1957,1959] was able to show that sonar 
amplitude fluctuations resulting from temperature measure- 
ments could be correlated to sea state. Additionally, he 
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showed that the linear dependence of the variance of sound 
pressure level on range [Sheehy 1950] was apparent- only when 
thermal microstructure was present. 

A modified geometrical optics technique has recently 
been proposed by Tatarski [1967]. The method is similar to 
the classical one except the path length is described as 
the eikonal or total path length, 6. However, one would 
assume that this would modify the restriction on the results 
to 



X « l 2 /0 (11) 

o 

No application of this method to the ocean has yet been made. 
2. Born Approximation 

The Born approximation is a perturbation technique 
applied to the wave equation which is valid only as long 
as the acoustic amplitude fluctuations remain small. This 
is not a serious restriction in the ocean, even under the 
most severe scattering conditions, but the model only accounts 
for a single scattering from a turbulent blob over the 
acoustic path length, L. This requires that 

X << a /aL (12) 

The geometrical optics results are a special, high frequency 
case included in this model. 
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Stone and Mintzer [1962, 1965 ] conducted a series 

of laboratory experiments from which they demonstrated 

solutions to the Born Model which predict the variance of 

2 

normalized sound pressure level, V , in a low frequency 
region (wave region) and in the ray theory region. For 
the wave region they showed: 

V 2 = | k Q 2 a 2 a L (Gaussian R ( p ) (13) 

V 2 = a 2 k 2 a L (Exponential R (p) (1*0 

o X 

v 2 = E[(^) 2 ] (15) 

o 

X << 2ua (16) 

X << 2ttL (17) 

k Q = 2ir/X (18) 

Where P Q = mean sound pressure level and AP = pressure level 
difference from the mean. For the ray region, no solution 
is possible if the exponential form of the spatial correlation 
function is used since the solution is divergent. But they 
showed that the following equation results from use of the 
Gaussian correlation function: 
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(19) 



V 2 = 4- /ir a 2 (L/a) 3 



Neubert and Lumley [1970] used data from the above 



experiments to compare with atmospheric data based upon the 
same model. Their results showed that the turbulent scales 
for air and water are of the same order of magnitude for 
"similar" turbulent flows. They were also able to better 
define the limit on wave length specified In relation (6), 
such that 



1 

X <. a/5[a k (u/c Q ) ] 2 



( 20 ) 




( 21 ) 



R t 2 • ISW 11 



( 22 ) 



where : 



X = Taylor turbulence microscale 
g 



o k = Prandtl number for the medium. 



u = local turbulent velocity 



c Q = mean speed of sound In the medium. 



R 



t 



Richardson number 



19 



and 



n 



g 



= Kolmogorov turbulence microscale. 



In addition, they found a restriction on acoustic frequency 
was necessary, though it was usually less restrictive than 
(20). This restriction is 

1 1 

f 1 750[a k (R t ) 2 ] 2 /A g (sec" 1 ) (23) 

It follows from (6) and (20) that 

1 

l Q z 0{a/50[a k (u/c o )] 2 } (cm) (2H) 



which lends additional meaning to relations (6), (7) and 

( 11 ). 

Whitmarsh, Skudrzyk and Urick [1957] conducted 
experiments in the ocean which largely agree with Stone’s 
and Mintzer’s laboratory results. However, they noted the 
variance in sound pressure level due to thermal scattering 
varied directly as the acoustic frequency and as the range 
squared in the Kolmogorov turbulence range. This result is 
obviously in the intermediate range between the ray and 
wave solution regions suggested by Stone and Mintzer. 

They also observed that the turbulent patches were 
not spherical, but were spheroidal in shape. In this case, 
a simple transform to cylindrical coordinates shows that the 
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Taylor microscale, a, must be replaced such that 

a 2 = (b 2 /3)(2 + c 2 /b 2 ) (25) 

where b = narrower dimension and c = longer dimension of 

the turbulent patch. This has the effect of reducing the 

focusing length of an individual turbulent "lens" which, 

in turn, increases the effective ray solution region. 

The following expression for the variance of the 

2 

sound pressure level, V , has been presented as valid over 
the entire wave and ray regions [Krasil’nikov and Obhukov, 
1956[: 

V 2 = k Q 2 a 2 a L[1 - l/Dtan _1 D] (26) 

D = TrL/k o a. (27) 

For this result, the wave region is defined when 

X » 2a 2 /L (28) 

and the ray region is defined for 

X << 2a 2 /L. (29) 
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3. Rytov's Method 

Rytov's method (method of smooth perturbations) is 
the method used by Tatarski [1961] to develop a stochastic 
wave model for the atmosphere. It includes the geometrical 
optics and Born approximation methods as special cases and 
includes multiple scattering effects as well. Unfortunately, 
the scattering effects are included such that they cannot 
be parameterized directly in terms of the turbulence. 

Either the spatial correlation or structure function 
may be used with this method as dictated by considerations 
regarding the stationarity of measurements to be analyzed. 

The restriction 



is necessary for reasonable solutions of the resultant 
equations. It is also necessary to compute the range of 
validity of the perturbation expansion when this method is 
used. This calculation will not be discussed, but Strohbehn 
[1968] suggested that the method is valid as long as the 
variance in the normalized log-amplitude of the sound pressure 
level is less than 0.8. He also argued that the results 
from sound phase calculations should be less restricted than 
those for amplitude. 

Black [1965] applied this model to certain data from 
the CHAIN experiment. One important conclusion he reached 




(30) 
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was that the standard deviation of the fluctuation in sound 



intensity through a turbulent medium, o a , and the principal 
period, T_ , of its fluctuation can be predicted by the equations 

S 



11 _ _7_ 

°s = 12 * 7 C n l12 x 12 • (31) 

T s = 7r/XT/2v n (32) 

C 2 n = C t 2 /4T 0 2 + C s 2 /1000S 2 + C y 2 /c q 2 = C t 2 /HT q (33) 

2 

C T 2 = D t (L)/L 3 (33a) 

2 

c s 2 = D S (L)/L 3 (33b) 

_2 

C y 2 = D v (L)/L 3 (33c) 



where : 

v n = absolute value of the medium's mean velocity 
across the acoustic path, 

T q = mean absolute temperature (°C), 

\ 

S = mean salinity in parts per thousand, 

C Q = mean speed sound (m/sec). 

Dm (L ) = value of the temperature structure function 
at range L, 

D g (L) = value of the salinity structure function 
at range L, 

D (L) = value of the turbulence structure function 
at range L. 



values for C n which 


Black 


observed ' 


In thermocline: 


5.2 


X 


io-V 1 


Near surface : 


1.0 


X 


io-V 1 


Below thermocline: 


2.1 


X 


io-V 1 
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4 . Debye Approximation 



The Debye approximation to the wave equation is the 
most general model currently available for the description 
of the turbulent ocean effects on sound transmission [Neubert 
1970]. It contains Rytov's method as part of its solution 
and it reduces to the Born model when the single scattering 
condition 

X >> 2 7T a/aZ ( 34 ) 

is applied. The resultant equation from the model, 

V 2 » 2 a 2 k Q 2 L e L (35) 

where L e = the Eulerian integral scale (e.g. the largest 
turbulent scale), has no special long range nor high frequency 
restrictions. In other words, it is good when the cumulative 
effects of small, continuous fluctuations of the local 
turbulent velocity strongly affect phase at longer ranges. 

No use has been made of this model in the ocean to date. 

C. MODEL APPLICATION 

Stewart and Grant [1962] conducted a turbulence experiment 
near the surface in the presence of large waves. They showed 
the Taylor microscale, a, was on the order of 6 cm which 
corresponds to the results (6.4 cm) Neubert and Lumley [1970] 
calculated from Sheehy's [1950] data. Since these measure- 
ments correspond well and fit the conditions of the current 
NPS experiment, the Neubert and Lumley values were used to 
provide a rough basis for further analysis. These values, 
coupled viith experimental parameters, are shown in Table 1. 
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The limits of applicability of each of the wave models and 

p 

the associated value of V calculated on the basis of this 

data is shown in Table 2. It would appear from the latter. 

table that the Debye approximation and perhaps the Born 

model for the wave region (if the limit imposed by relation 

(6) is not too strong) would be most applicable in this 

2 

case. The large, difference in V resultant from the Born 
model compared to that from the Debye model could be the 
result of poor definition of the micro and integral scales. 

One should also note that the acoustic near-field 
criteria for a piston source imposes an additional requirement 
on the wavelength: 

A > r 2 /L (36) 

where r = radius of the transmitter active disc. Since 
relation ( 36 ) requires A > 0.6 cm, geometrical optics and 
Rytov’s method can be immediately discarded. 
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TABLE 1 



PARAMETER VALUES FOR MODEL TESTING 



’arameter 


Assumed Value 


Source 


°o 


1505 m/sec 


Experiment 


u/0 o 


7 x 10" 5 


NL 


L e 


40m 


NL 


\s 


6.4 cm 


NL 


R t 


6300 


NL 


CT k 


7 


NL 


2 

a a 


5.0 x 10 -7 


S 


L 


200 cm 


Experiment 


f 


20 khz 


Experiment 


A 


7.5 cm 


Experiment 


C n 


-4 -1 
1.0 x 10 m x 


B 


u 


10.5 cm/sec 


C 


n g 


258 cm 


Equation 22 


a 


2.4 cm 


Equation 21 


*0 


2 . 2 cm 


Equation 24 


a 


4.6 x 10 - * 1 


C 


k o 


0 . 84 cm - ' 1 ' 


Equation 18 


D 


130 


Equation 27 


r 


10 . 88 cm 


Experiment 



NL = Neubert and Lumley [1970] 

S = Sheehy [1950] 

B = Black [1965] 

C = Calculated from above data. 
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TABLE 2 



RESULTS FROM MODEL TESTING 



Model 


£ 




Most Restrictive Limit 


Geometric 


5.8 x 10 -2 


(19) 


_2 

A « 1.1 x 10 cm 


(7) 


Bom 

(Gaussian fom) 


6.3 x IQ' 5 


(13) 


-2 

1.0 x 10 « A « 


15cm(12,l6) 


Bom 

(Exponential form) 


7.1 x 10" 5 


(14) 


_2 

1.0 x 10 « A « 


15cm(12,l6) 


Bom 

(Neubert Limits) 






A £ 20.9cm (18) 




Bom (K-O) 


7.1 x 10" 5 


(26) 


A » 2.9 x 10 _2 cm 


(28) 


Rytov 


1.1 x 10" 2 


(3D 


A « 0.5 cm (30) 




Debye 
k = 40m) 


2.4 x 10" 1 


(35) 


None 




Debye 
(L e = 3 m) 


1.8 x 10 -2 


(35) 


None 





HOMOGENEOUS ASSUMPTION 
A << 2 .2 cm ( 6 ) 

ACOUSTIC FAR- FIELD CRITERION 
A > 0.6 cm (36) 

Note: Numbers in parenthesis refer to applicable equation. 
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III. EXPERIMENT AND DATA ANALYSIS 



A. EXPERIMENT 

Measurements were taken from May 24 to 26, 1974 from 
the Naval Undersea Research and Development Center’s (NUC) 
oceanographic research tower located 1.6 km off Mission 
Beach San Diego, California. The tower is a stable, three- 
tiered concrete and steel structure built into sandy bottom 
in water of 18.5 m mean depth. It provides enclosed spaces 
for electronic equipment and open-deck working areas for 
handling instruments. Winch controlled trolleys are avail- 
able on three sides of the structure on which devices may 
be lowered to within 1 m of the bottom. Instruments for 
this experiment were installed on an aluminum and steel 
frame attached to a trolley on the west side of the tower 
which faced the predominant swell. A schematic view of 
the experiment layout is shown in Figure 1. 

Basic meteorological and oceanographic data were avail- 
able from equipment permanently installed on the tower. A 
bathythermograph trace obtainable from the north side of 
the structure facilitated placement of the instrument frame 
in the desired vertical position with respect to the 
thermocline . 

The general environmental conditions encountered were 
moderate winds, sea and swell from the West. Light winds 
to 5 kts, overcast skies and swell to 1 m were usual in the 
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morning (see Table 3). In the afternoon, winds increased 
to 15 kts and seas of 0.75 m with accompanying swell of 1 m 
was common. The diurnal variation in thermocline depth was 
approximately 12 m. 

1. Instruments 

Temperature measurements were made with Penwell 
Electronics Inc. type k*J96, double-bead thermistors which 
exhibit low self-heating and well-matched, high-resolution 
characteristics. Calibration of the thermistors, together 
with associated circuitry, was conducted in the laboratory 
under closely controlled, stable conditions. A Hewlett- 
Packard HP-2801A quartz thermometer with an accuracy of 
± .006°C was used as a temperature standard and a decade 
box was used as a resistance standard. The same decade box 
was used to calibrate the thermistor array in the field to 
eliminate resolution difficulties resulting from the poor 
overall accuracy of the decade box. 

Water particle velocities were measured using Marsh- 
McBirney electromagnetic flowmeters (models 721 and 722). 
Since these flowmeters were calibrated by the manufacturer 
only for steady flow conditions, each probe was recalibrated 
for unsteady flow by oscillating it in a water tank at rates 
corresponding to expected experiment velocities. 

A Baylor Co. wave staff system, model 13528R, was 
used to measure sea surface elevation. Calibration was 
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26 Apr 13.6 13.4 12.65 11.6 10.9 10.7 10.7 



accomplished at the NUC tower by successively shorting out 
sections of the sensor active element while measuring output 
voltage . 

A Naval Research Laboratory, Underwater Sound 
Reference Division (USRD) type F-33 general purpose uni- 
directional transducer was used as a CW sound transmitter. 

It features two piezoelectric ceramic arrays mounted coaxially 
which, when excited in parallel, exhibit exceptional direc- 
tivity. This was desired to reduce interference from surface, 
bottom and equipment frame reflections. The beam widths 
(-3db down) for 20kHz and 40kHz were 20° and 10° respectively. 
USRD type F-50 lead zirconate-titanate omnidirectional 
transducers were used as receivers. 

The sonar array was calibrated in an anechoic test 
tank filled with fresh water and the hydrophone calibrations 
compared satisfactorily with calibrations previously per- 
formed at USRD. . The frequency standard used was a Hewlett- 
Packard HP-2801A quartz thermometer configured as a frequency 
meter and the voltage standards were an RMS calibrated Fluke 
digital voltmeter and a Hewlett-Packard model HP 3400A RMS 
voltmeter. Phase measurements were conducted using a Dranetz 
model 305 phase meter. Calibration and accuracy measurements 
for the peripheral equipment associated with the sound arrays 
were taken concurrently with transducer calibration. 

The instrumentation accuracy, response and limitations 
are summarized in Table 4. 
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INSTRUMENT SUMMARY 
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2 . Instrument Arrangement 



Two arrays of thermistors were mounted on one side 
of the equipment frame in a common vertical plane. Seven 
thermistors were mounted along a horizontal member and 7 
were mounted along a vertical member as shown in Figure 2 . 

The instrumentation and abbreviated notation used in the 
text is given in Table 5. One thermistor in the vertical 
array was common to the horizontal array. The thermistor 
probes were spaced in an approximate geometric progression 
so that a maximum number of spatial combinations could be 
achieved for adequate definition of the spatial correlation 
and structure functions with a limited number of probes. 

The minimum spacing ( 5.0 cm) was determined by the physical 
size of the thermistor mounts. This was compatible with 
the smallest wavelengths which could be measured within the 
time constant of the thermistor. Maximum thermistor spacing 
was determined by physical constraints on the length of the 
support pipe to prevent strumming. 

Three velocimeters were placed so that the horizontal 
and vertical particle velocities could be measured in the 
same plane as the thermistor arrays. Additionally, they 
were placed in a vertical line directly under the Baylor 
wave gage with the horizontal component measured from the 
west. The temperature and particle velocity devices were 
offset 29 cm from the support frame by mounting brackets to 
minimize false turbulence effects generated by the frame. 
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FIGURE 2. RELATIVE POSITIONS OF INSTRUMENTS ON FRAME (SIDE VIEW) 
INSERTS SHOW SONAR ARRAY CENTERLINE POSITION RELATIVE TO FRAME BAYLOR 
WAVE GAGE WAS LOCATED DIRECTLY ABOVE VERTICAL ARRAY DEVICES. 




TABLE 5 

CHANNEL ABBREVIATIONS 



Data Channel 
Vertical Temperatures 
Horizontal Temperatures 
Sound Amplitude 
Sound Phase 
Sea Surface Elevation 
Horizontal Plow Velocity 
Vertical Flow Velocity 
Conductivity 



Abbreviation 

TMV1-TMV7 

TMH1-TMH7 

SAM1-SAM3 

SPH1-SPH2 

WAVE 

FVAY-FVCY 

FVAX-FVCX 

SAL1-SAL2 



Order in Fig. 2 



From Bottom 

From West 

From East 
or 

_From Top 



From Bottom 



From Bottom 
From West 
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The sonar array was separately mounted on a tripod 
base which pivoted vertically around the main equipment 
frame but which was held rigid by clamps during recording 
periods. This arrangement facilitated horizontal and 
vertical sound path measurement (Figure 2). The centerline 
of the transmitted sound beam was 57 cm away from the 
vertical plane of the thermistor arrays to minimize unwanted 
reflections from the other instruments. 

Three F-50 receivers were placed along the center 
of the transmitted beam outside the near field at 2.0, 2.3 
and 2.9 ni from the transmitter face. The minimum separation 
(30 cm) was selected to maximize the number of wavelengths 
through which the lowest frequency (20kHz) sound would 
travel between receivers and to remain near the lower spatial 
scales of the measured temperature. The maximum sound path 
length of 2.9 m was selected to minimize reflections from 
the frame and mountings . 

3. Pre-recording Signal Processing 

The signals from the wave height and particle velocity 
devices were recorded directly. All other signals were 
processed through a combination of filters, amplifiers and 
dc bias removal devices to maximize signal-to-noise (S/N) 
ratio before recording. 

a. Temperature Circuitry 

Each of the horizontal and vertical temperature 
arrays was connected to a separate temperature array drawer 
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which contained Wheatstone bridge networks, amplifiers and 
calibration facilities for seven channels (see Figure 3). 

The Wheatstone bridge network provided the capability to 
match the response of each channel. Additionally, balance 
adjustments were provided to adjust the operating point of 
the amplifiers to the mean temperature so that distortion 
during recording was minimized, 
b. Sonar Circuitry 

As shown in the block diagram (Figure *0 , signals 
from the F-50 transducers were amplified and then band-pass 
filtered at the operating frequency. The signal was then 
processed through a second-detector-type demodulator. The 
dc bias was removed from the detected signal and the ampli- 
tude modulation (AM) signal was then amplified and low-pass 
filtered. 

Dranetz type 305 phase meters were connected 
between the band-pass filter outputs from receiver hydrophones 
one and two and, similarly, between hydrophones two and three. 
The dc bias was removed from the output of each phase meter 
and the signal was then amplified and low-pass filtered. 

Before each run the dc bias setting was made and 
the rms voltage from each band-pass filter output and mean 
phase shift read from each phase meter was recorded. This 
technique allowed much higher signal amplification to be 
used prior to recording because the signal amplifiers were 
not saturated by the large dc segment of the signal. 
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OCEAN ENVIRONMENT 



r 




FIGURE 3. BASIC TEMPERATURE MEASUREMENT 

CIRCUIT BLOCK DIAGRAM. ONLY I OF THE 7 CHANNELS IS 
SHOWN IN THE TEMPERATURE ARRAY DRAWER FOR CLARITY. 
TWO IDENTICAL ARRAY DRAWERS WERE USED DURING THE 
EXPERIMENT. 
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SONAR ARRAY IN OCEAN ENVIRONMENT 
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FIGURE 4. BASIC SONAR PHASE AND AMPLITUDE 

MEASUREMENT BLOCK DIAGRAM. TWO PHASE CHANNELS AND 
THREE AMPLITUDE CHANNELS WERE USED, HOWEVER ONLY 
ONE OF EACH IS SHOWN FOR CLARITY. 
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Experimental Procedure 



Twelve data recording runs of 29 to 73 minutes 
duration were taken at various depths above, in and below 
the thermocline. A total of 28 channels of data were 
recorded for each run including 2 channels dedicated to 
experimental conductivity devices [Kane 197*1]. Surface 
elevation, vertical and horizontal particle velocities, 
temperature, sound amplitude modulation, sound phase modu- 
lation and conductivity were measured simultaneously. 

Run depths were selected on the basis of a BT trace 
taken Just before the run, but the final run depth was 
chosen to optimize the degree of fluctuation in temperatures 
monitored from the vertical thermistor array. The sonar 
array was held in a fixed position during each recording 
run. Runs were made with the sonar array oriented hori- 
zontally, vertically and at an angle with respect to the 
support structure. This allowed data to be obtained at 
various attitudes through and along the thermocline. A 
summary of all twelve runs is given in Table 6. 

A Sabre III and a Sangamo model 3500, both 1*1 channel 
FM tape recorders, were used to record the analog data. A 
timing pulse was recorded on one channel of each recorder 
at the beginning of each run to provide temporal continuity 
between recorders. 

The analog recordings were digitized using a hybrid 
system including a Scientific Data System model 5000 analog 



SUMMARY OF RUNS 
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computer and a Xerox Data Systems model 9300 digital computer. 
Digital data preparation and analysis was performed on an 
IBM model 360 computer. The digitizing data rate used was 
25Hz, however 60Hz electronic noise was coupled into the 
signal, so low-pass filters were used to remove it. Addi- 
tionally, noise spikes which existed in some channels were 
removed by a digital editing routine. Data handling capacity 
was improved by reducing the number of samples by five using 
an inverse transform, numerical low-pass filter [Davidson 
1970]. The latter reduced the analysis Nyquist frequency 
to 2.5Hz. The linear trend and remaining dc bias were also 
removed before analysis. 

Spectral analysis was performed on 26.6- minute 
segments of certain data sets. This choice dictated a 
minimum detectable frequency of 0. 00625Hz. A standard lag 
time of 5% of the record length was chosen which resulted in 
40 degrees of freedom for each computed spectrum. The 
resultant 80 % confidence limits of the chi-squared distri- 
bution were between 0.73 and 1.30 of the measured spectral 
estimates. A Parzen window was applied to the correlation 
function during the spectra calculations to minimize the 
effect of finite record length on the spectral analysis. 
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IV. EXPERIMENT RESULTS 



During this experiment, the primary factors influencing 
the results were the sea state, the thermocline gradient and 
the position of the instruments with respect to the thermo- 
cline. The sea state remained nearly constant throughout 
the recording periods although the peak of the swell showed 
a gradual shift to higher frequencies. The thermocline 
depth and gradient continuously fluctuated under the influence 
of internal waves which are characteristic of the waters off 
San Diego during this time of year. In order to resolve 
some of the effects of the thermolcine variability, analyses 
in and above the layer were required. Experiment runs one 
and two were chosen for this purpose. Additionally, run 
12 was examined to check the results since: 1) it included 

data in and above the thermocline in a single record; 2) the 
thermocline gradient was intermediate between the gradient 
encountered in runs one and two; 3) the sonar array was 
angled toward the vertical rather than horizontal during the 
run and 4) the run was conducted two days after runs one 
and two under slightly different sea conditions. 

A. ENVIRONMENT 

The thermocline migrated diurnally from near the surface 
in the morning to about l^m in the evening. It was generally 
diffuse near the surface but near the bottom it was compressed 
so that a very sharp thermal gradient (2°C/m) was created. 



Superimposed upon this effect was an apparent four-hourly 
dilation and contraction of the isotherms. 

Figures 5-9 show representative bathythermograph (BT) 
traces taken periodically during the experiment. The traces 
are plotted so that the 12°C mark corresponds to the time 
the BT was taken. Recording run envelopes in depth and time 
are included to show the position of the instrument array 
relative to the thermocline. The 11° , 12° and 13°C isotherms 
are shown to indicate the relative intensity of the tempera- 
ture gradient during the runs. Figure 9 is an expanded 
time plot for BT's taken during run 12 which shows a large 
internal wave having a 36 min period; this was observed 
throughout the experiment. 

The four-hourly variation in the thermocline gradient 
probably was caused by internal wave interference effects 
since analysis showed that these variations were associated 
with periods of high temperature variance but not with sur- 
face waves nor heat budget considerations. 

Surface slicks presumably associated with internal waves 
were seen propagating from the West at certain intervals 
during the experiment. A typical example of the effect of 
these waves is shown in the vertical temperature analog 
records of run number 12 (Fig. 10). The large 36 min wave 
peak passing the vertical temperature array is shown by the 
dramatic temperature decrease which occurred at about 1800. 
Nine and 1.5 min period, small-amplitude internal waves 
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FIGURE o. BATHYTHERMOGRAPH TRACES FOR 24 APRIL 1974. TYPE OF THE TRACE IS LOCATED 
AT THE I2°C TEMPERATURE MARK. RECORDING RUN ENVELOPES IN DEPTH AND TIME ARE INCLUDED 
TO SHOW THE POSITION OF THE INSTRUMENT ARRAY RELATIVE TO THE THERMOCLINE . THE 11°, I2°AND 
13° C ISOTHERMS ARE SHOWN TO SHOW THE RELATIVE INTENSITY OF THE TEMPERATURE GRADIENT 
DURING THE RUNS. 
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FIGURE 6. BATHYTHERMOGRAPH TRACES FOR THE MORNING OF 25 APRIL. 
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FIGURE 7 . BATHYTHERMOGRAPH TRACES FOR THE AFTERNOON OF 25 APRIL. 
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FIGURE 8. BATHYTHERMOGRAPH TRACES FOR 26 APRIL. 
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FIGURE 9 . SHORT TIME INTERVAL BATHYTHERMOGRAPH TRACES FOR THE 
PERIOD OF RECORDING RUN 12, 26 APRIL. 



effects are seen superimposed upon the larger internal wave. 
The "grass” is surface swell related. Spectral analysis 
showed that a broad band of internal wave periods from 
160 secs (maximum detectable by the method) to as low as 
30 secs were superimposed on the large wave as well. 

Figure 11 shows the resultant sound effects from the 
internal wave temperature perturbations of Fig. 10 during 
run 12. The sound array was angled 25° to the vertical in 

run 12. The sound amplitude modulated signals SAMI and SAM2 

were initially above the layer and were not affected by 
the turbulent thermocline losses as severely as SAM3 (below 
the layer). As the layer moved up, the interference and 
diffraction effects of the turbulent blobs [Whitmarsh and 
others 1957] in the vicinity of SAM3 were diminished. This 
is shown by the increase in sound pressure level (about 
10 J ub) on the SAM3 analog record. 

When the layer moved into the vicinity of SAMI and SAM2, 

the sound pressure level at these transducers was reduced 

by about 6 x 10 yb and 3 x 10 yb , respectively. The 
thickness of the "grass" on these plots is an indication of 
the wave-induced amplitude scintillation effects (~2 to 5 
x 10 yb) caused by pumping the turbulent blobs back and 
forth across the sound path. (Similar results occurred when 
the sound path was horizontal, but the scintillations were 
then caused by waves pumping the thermocline up and down. ) 
Additionally, the movement of the layer (AT = .3°C) between 
SAMI and SAM2 caused sound phase modulation, SPH1, to phase 
shift about 30°. 
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B. TEMPORAL STOCHASTIC ANALYSIS 

Time series analyses were performed on 26.67 minute 
segments of each of the sound channels, the wave record, 
three temperature records and two water particle velocity 
records. Analyses included for each signal the calculation 
of the first and second moments, variance, mean period 
(calculated from total zero upcrossings) , autocorrelation 
function and autospectrum. Also, the crosscorrelation, 
quadrature and co-spectra, phase difference, coherence- 
squared and linear regression analysis were calculated 
between specific pairs of channels (see Table 7). Examples 
of cross-spectral analysis results for runs two and 12 are 
shown in Figs. 12-15. Additional results are given in 
Appendix A. Thermistor 2 on the horizontal array, TMH2 , 
was selected for temperature comparisons with the sound 
records because this thermistor was most central to the 
sound propagation path during the runs. 

1. Basic Statistics 

Representative basic statistics of temperature, 
waves, sound amplitude and sound phase are shown in Table 
8. By examining the order of magnitude of these statistics 
in relation to the depth of the run and the proximity to 
the thermocline, certain qualitative relationships between 
the parameters can be seen: 

1) The magnitude of the temperature variance was 
determined primarily by proximity to the layer and its mean 



Ordinates For Regression Analyses 



TABLE 7 

CHANNELS COMPARED IN TEMPORAL STOCHASTIC ANALYSES 



Channels 

SAMI 

SAM2 

SAM3 

SPH1 

SPH2 

TMH2 



Abscissas For Regression Analyses 



TMH2 WAVE SAL2 

XXX 

X 

X 

XXX 

X 

X X 
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PFTflSF COHERENCE LOG SPECTRRL DENSITV 

OOO +180 0 1-3 -5 -3 






FIGURE 12. Run 2 TMH2 vs WAVE Spectra 
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rHHbt'. COHERENCE LOG SPECTRRL DENSITY 
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. Run 2 SAMI vs TMH2 Spectra 
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FIGURE lk 
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* Anomalous because sonar array angled 25° to vertical vice horizontal. 



period was controlled by surface wave swell until the layer 
was very close to the thermocline. 

2) Sound amplitude variance was controlled by proximity 
to the thermocline (temperature gradient) when the sonar 
array was horizontal. The movement of the thermocline 

with respect to the sound path must therefore be the primary 
influence on sound amplitudes in this configuration. 

3) The amplitude variance when the array was nearly 
vertical was much smaller than when it was horizontal. 

This indicated the effect of turbulence proximity to the 
hydrophone was smaller than the effect of temperature 
gradient motion over the sound path. 

4) Since the mean period of the sound fluctuations was 
much smaller than the mean period of the swell at depth, 
turbulent interactions related to the higher frequency 
internal wave motion of the thermocline were emphasized 

in the sound measurements. 

5) When the sound array was vertically oriented vice 
horizontal, much larger phase variations occurred even 
though amplitude variations were smaller. This resulted 
from the large temperature variations over thepath length 
caused by waves pumping the thermal gradient up and down. 

This also showed that the temperature variation effects 
on sound amplitude and phase were in quadrature. 

6) The order of the mean sound phase variances indicates 
that phase was primarily influenced by surface waves when the 
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array was horizontal. Since spectral analysis showed that 
the wave coherence with temperature was very low at the 
surface wave frequencies, this effect may be bubble related. 

2 . Spectral Analysis 

Cross-spectral analysis of the temperature, waves 
and sound records for the horizontal sonar runs one and 
two showed that the sound amplitude changes (Figs. 13 and 
14) were in phase with the waves and l80° out of phase with 
temperature in the regions of highest coherence. This was 
expected in the stable density structure encountered. 
Additionally, the coherence squared function (simply re- 
ferred to as coherence) showed that only the long-period, 
prominent swell was important to the sound amplitude modu- 
lation because the higher frequency wave particle velocities 
were much more rapidly attenuated with depth. 

Sound phase (Fig. 15) was found to vary directly 
with temperature at the internal wave and surface wave 
frequencies although coherence at the surface wave frequen- 
cies was much smaller than at the internal wave frequencies. 
This indicated that the turbulence associated with the layer 
was more important than bubbles or surface wave perturbation 
of the thermocline in influencing sound speed over the 
horizontal sound path. A comparison of the coherences of 
Figs. 12 and 13 shows that the temperature effect is at 
least as important as bubbles even at the surface wave 
frequencies . 
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Similar cross-spectral analyses with the sonar path 
angles 25° from the vertical (run 12) yielded markedly 
different results. For the sound path entirely above the 
layer (SAMI), coherence with temperature was small over all 
frequencies. Small coherence peaks occurred at 0.15Hz and 
0.7Hz which indicated that the vertical motion of the turbu- 
lent blobs forced by surface swell had some effect in the 
angled sound path. 

Sound phase coherence with temperature was largest 
at the internal wave frequencies but smaller, broad-banded 
coherence at the surface wave frequencies was also observed. 
This confirmed that the phase was being controlled by the 
vertical motion of the thermocline. 

Most of the sound path associated with SAM2 (2.3m) 
was not in the thermocline during the majority of run 12. 
However the hydrophone was near the top of the layer. As a 
result, its amplitude was in phase with temperature at the 
internal wave frequencies and coherence at other frequencies 
was almost zero. This verified that the interference effect 
of the proximity of the thermocline-associated turbulence 
to the receiver was theprimary effect on SAM2 variability. 

Very high coherence and 180° phase difference between 
TMH2 and SAM3 at the internal wave and the surface wave 
frequencies was observed. Since the acoustic path associated 
with SAM3 generally spanned the thermocline, this indicated 
that the variability of amplitude across the entire 
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thermocline was controlled by refraction and diffraction 
similar to the horizontal case. The effect of diffraction 
was smaller than for the horizontal case, because the len- 
ticular shape of the turbulent blobs present broader profiles 
to the vertical sound path. The reduced diffraction effects 
caused the variance in sound amplitude for run 12 to be an 
order-of-magnitude less than that for run 1. 

3. Regression Analysis 

2 

Sound amplitude and phase variances, V (SAM) and 

2 2 

V (SPH), were plotted versus TMH2 variances, V (T) , and 

least squares regression analyses were applied to the data. 

2 2 

Linear relationships between V (SAM) and V (SPH) versus 
2 

V (T) were found to exist for runs one and two but not for 
run 12. This was because of insufficient variability in 
the portion of run 12 analyzed. The empirical relationship 
below was developed from the mean regressions calculated 
for runs one and two: 

V 2 (SAM) = (0. 34yb 2 /°C 2 )V 2 (T)e _L/x (37) 

where L = acoustic path length and x = the distance required 
to reduce V (SAM) by 1/e (determined empirically). x is 
near the size of the Taylor microscale. 

The derivation of Equation (37), plots of the 
variances over the run segments and the regression plots 
are given in Appendix B. 
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C. SPATIAL ANALYSIS 

Temperature spatial correlation and structure functions 
were calculated directly from ensemble averages of 26.67 min 
records from all possible combinations of thermistor pairs 
along the vertical thermistor array. Twenty-two spatial 
lags with values between 0 and 260 cm were thus determined 
from seven thermistors. The value of (see Equation 33a) 
and the Taylor microscale as well as the Eulerian integral 
scale were derived from the above functions (Table 9)» 

However, since large variations were seen in these functions 
the validity of the calculated scales was questionable. 

In an attempt to remove the non-stationary effects 
resulting from internal waves, analyses using ensemble 
averages of short segments having nearly constant turbulent 
intensity were conducted for runs two and twelve. The 
segments were chosen so that examples of low, medium and high 
turbulent intensity could be examined (Table 10). Because 
the thermocline steepness was highly variable along the 
vertical thermistor array, however, only the medium inten- 
sity correlation function was smooth enough to provide satis- 
factory Eulerian integral scales and the structure function 

of all samples varied widely. Values of the Taylor microscale, 

2 

Eulerian integral scale and generated from hand-smoothed 

results are given in Table 11. The correlation function and 
structure function values as well as the mean variance for 
the vertical array thermistors are given in Appendix C. 
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TABLE 9 



Temperature Spatial Scales (cm) for 26.67 min Ensembles 



Run if 


Taylor 1 

Microscale 


Taylor 2 

Microscale 


Eulerian 
Integral Scale 


C\J 

cF 


1 


15 


40 


36 


1.45x10"' 


2 


178 


158 


175 


3.6xlO“ 3 


12 


60 


100 


36 


8.3xlO“ 3 



TABLE 10 

Relative Turbulent Intensity for Special Spatial Analyses 



Qualitative 

Intensity 


Ensemble 
Size (min) 


TMH2 Mean 
Variance (°C) ^ 


Range 

(°C)2 


Run 


Low 


4.6 7 


0.27xl0 _1 


±0 . 77xl0 -2 


12 


Medium 


4.67 


4.24xl0 _1 


±4.95xl0 _1 


2 


High 


3.13 


7.79xl0 _1 


±2. 2 3x1 0" 1 


2 



Temperature 

Qualitative 

Intensity 


TABLE 11 

Spatial Scales (cm) for Special Ensembles 
i 2 

Taylor Taylor Eulerian ^ 2 

Microscale Microscale Integral Scale °T 


Low 


23 


6.8 


24 


3.3x10"- 


Medium 


71 


40 


35 


4 . 4x10" - 


High 


15 


20 


63 


5.4x10"- 



Taylor 1 = Taylor microscale from spatial correlation 

p 

Taylor = Taylor microscale from structure function 
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D. APPLICATION OF MODELS 

While it is apparent that some of the scale estimates 
calculated from the vertical thermistor array may be erron- 
eous, it is instructive to apply the theoretical models 
for the sound amplitude variance to see which model is 
most applicable in the non-stat ionary case. Table 12 lists 
the equations necessary to calculate the basic parameters 
for application of the models. For the purposes of this 
examination, the values of the Richardson number (R^)> the 
Prandtl number (a ) and u/c were assumed to be the same 

K O 

as reported by Neubert [1970] (see Table 13 for values). 

The speed of sound is used to calculated Wilson’s Equation 
[Urick 19671. 

Tables 13 and 14 list the resultant parameters of each 
sample for which scales were determined. It is noteworthy 
that the values of the refractive index (a) calculated above 
are an order of magnitude above the near surface results 
reported by other investigators [Haley 1972, Medwin 1973] 

(see Table 15) except in the low turbulent intensity cases. 

This is further evidence that the relative influence of the 
turbulent thermocline in sound scattering is much greater 
than bubbles. 

If the order of magnitude of the Taylor microscale is 

correct (see Table 13 for comparative values) the values of 
2 

a a are one to four orders of magnitude greater than px-eviously 
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TABLE 12 

EQUATIONS USED FOR MODEL PARAMETER EVALUATION 



Parameter 



c 



o 



a 



C n 



Equation 

c Q = 1449 + 4 . 6T - 0.55T 2 + 0.17d 
a 2 = [4. 6 2 V 2 (T) + 0.055 2 V 4 (T)]/c o 2 
X = a /a, = 2.65a 

g K 

n g = A g 4 (H t /15) = 2.07 x 10 V 

C n 2 = C T 2 / 4 T o 2 

l Q * = 0(2. 2/2. 4a) = 0(0. 9a) 



T = mean temperature (°C) 

d = depth (m) 

2 2 
V (T) = variance in temperature (°C ) 

a = Taylor microscale (cm) 

n = Kolmogorov microscale (cm) 

o 

* equation for l assumes u/c is a constant 

o o 
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TABLE 13 



Experiment Parameter Values for Model Tests 



Parameter 


RUN 1 


RUN 2 


RUN 12 


c Q (m/sec) 


1516 


1498 


1515 


a* 


4 . 5xl0 - ^ 


2 .lxlO -3 


7.9x10 


L ( cm) 


200 


200 


200 


a ( cm) 


15 


178 


60 


V 


13.5 


160 


54 


X ( cm) 


7.5 


7.5 


7.5 


k Q (cm -1 ) 


0.84 


0.84 


0.84 


C n 


1. 45x10 -I< 


1 . 05xl0 - ^ 


2.4x10 


L e (cm) 


36 


175 


36 


°k 


7 


7 


7 


X (cm) 

o 


40 


472 


159 


n g (m> 


10.5 


21.4 


19.2 


R t 


6300 


6300 


6300 


a 2 a 


3.0xl0 -6 


-4 

7.9x10 


3.7x10 



-4 



-4 
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TABLE 14 



Short Sample Parameter Values for Model Tests 



Parameter 


Low 

Intensity 


Me dlum 
Intensity 


High 

Intensity 


c Q (m/sec) 


1512 


1508 


1500 


a 


5 . OxlO - ^ 


2. 0xl0 -3 


2.7xl0 -3 


L ( cm) 


200 


200 


200 


a ( cm) 


23 


71 


15 


* 


21 


64 


13 


X (cm) 


7.5 


7.5 


7.5 


k Q (cm -1 ) 


0.84 


00 

0 


^r 

OO 

• 

O 


C n 


1.05x10"^ 


I.l6xl0 -4 


1 . 29xl0 -i * 


L e (cm) 


24 


35 


63 


°k 


7 


7 


7 


X g (cm) 


6l 


188 


40 


n g (cm) 


2.8 


8.5 


H 

• 

OO 


R t 


6300 


6300 


6300 


2 

a a (cm) 


5. 8xl0 -6 


2 . 8xl0 - ^ 


1 . lxlO - ^ 
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Comparative Values of a, a and aa 
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reported. This result has a direct relationship (in the 
Born model) to the variance of sound amplitude. 

The variances of sound amplitude calculated by each 
model previously discussed (Geometric, Born, Rytov, Debye) 
are presented in Tables 16 - 19. In general, the best fit 
to the measured sound variances were calculated from the 
Born model using the Krasil ’ nikov-Obhukov formulation. A 
surprisingly good fit to the true variance was also obtained 
from the geometric model in the mid-intensity range. Because 
of the above, the stimates of the Taylor microscale seem 

reasonable . 

2 

Since the accuracy of the scales, C,p and L e , the Eulerian 
integral scale, is in question, no firm statement can be 
made about the adequacy of the Rytov and Debye models. 

However, because L e is thought to be correct for the mid- 
intensity ensemble, the Debye model appears incorrect. 
Presumably, calculation of the structure function and corre- 
lation function from the thermistors in the horizontal array 
over stationary ensembles would resolve this question. 
Unfortunately, this was not accomplished in time for 
inclusion in this thesis. 
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TABLE 16 



Geometric Model Experimental Results 
V 2 = Vl5 /? ct 2 (L/a) 3 ; 7.5 << l 0 2 /L (cm) 



Ensemble 


Model V 2 


Observed V 2 (yb) 2 


\ 

OJ 

o 


Run 1 


8xl0 _i< 


1.12x10"^ 


0.91 


Run 2 


3xlO“ 6 


1.^5xl0“ 3 


128 


Run 12 


1. lxlO -5 


A.05x10“ 5 


1^6 


Low Intensity 


7. 8xl0 -5 


7.1xl0“ 6 


2.2 


Mid Intensity 


4.2xl0“ 5 


9.5xl0" 5 


21 


High Intensity 


8.2xl0 _3 


3. 0x10”^ 


0. 85 




TABLE 


17 




Born Model Experimental Results 
[Krasil'nikov and Obhukov] 




O ? ? 

V = k a aL(l 


- 1/Dtan 1 D) 


; D = TTL/k Q a 




Ensemble 


Model V 2 


p 

Observed V 


D 


Run 1 


. 3x10”^ 


1.12x10 


50 


Run 2 


7. 9x10”^ 


1.45xl0" 3 


OJ 

• 

-=T 


Run 12 


3. 7^xl0 -5 


.05xl0" 5 


12.5 


Low Intensity 


5. 8xl0“ 6 


7.1xl0 -6 


32.5 


Mid Intensity 


2 . 8x10”^ 


9. 5xl0" 5 


10.5 


High Intensity 


1 . 1x10”^ 


3 . 0x10”^ 


50 




D >> 1 


Wave Region 





(cm) 
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TABLE 18 



Rytov Model 

V 2 = 161 C n 2 L 11/6 -7/6 



Ensemble 


Model V 2 


Run 1 


5.2xl0 -3 


Run 2 


2. 8xl0~ 3 


Run 12 


1. 5xl0 -2 


Low Intensity 


2. 8xl0~ 3 


Mid Intensity 


3. 4xl0 -3 


High Intensity 


4.2xl0“ 3 



Experimental Results 



7.5 « V£ o */L 
2 

Observed V 


3 A 


1.12X10” 1 * 


5 


1 . 45x10” 3 


150 


4.05xl0" 5 


35 


7.1xl0 -6 


9 


9.5xl0" 5 

-4 




3.0x10 4 


5 



TABLE 19 

Debye Model Experimental Results 



2 2 2 
V ‘ 2 « k o L e L 



Ensemble 


Model V 2 


Observed 


Run 1 


2 . Ixl0~ 3 


1 . 12x10”^ 


Run 2 


0.22 


1.45x10 -3 


Run 12 


6. 3xl0 -3 


4.05xl0 -5 


Low Intensity 


1. 7xlO -3 


7.1xl0 -6 


Mid Intensity 


4xl0 -2 


9.5xl0" 5 


High Intensity 


0.13 


3.0xl0" 4 



74 



V. CONCLUSIONS 



The theories of turbulence which have been developed 
for the ocean require some degree of stationarity in the 
medium before they may be applied. This paper has examined 
four basic stochastic models of acoustic wave propagation 
under the most severe conditions of non-stationarity , e.g. 
near a sharp thermocline which is being oscillated by both 
large surface swell and by internal waves. 

The experimental conditions specified that the solution 
form would be in the wave region for which the Born model 
would be most accurate. The solution of the Krasil 'nikov- 
Obhukov form of the Born model from directly measured, 
ocean parameters gave reasonable results. 

The solution of the Rytov and Debye models for the same 

conditions yielded results two orders of magnitude higher 

than was actually measured. Because evaluation of scales 

was difficult, this was expected even though both models 

are more general than the Born model and, in fact, contain 

its solution in a subset. This also showed the determination 

of the Taylor microscale was more reliable than the deter- 

2 

mination of the structure constant, C,p , or the Eulerian 
integral scale for measurements across the thermocline. 

An attempt was made to apply the correlation function 
and structure function piecewise to small sections of data 
which were thought to be relatively stationary. For 
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application across the thermocline, this method did not work 
well because selection of such sections is very difficult. 

The results of the time series analyses showed that the 
losses over the acoustic path were greatest when the array 
was horizontal near the thermocline because the diffractive 
losses were an order of magnitude greater than those for 
the nearly vertical case. Interference effects due to the 
proximity of turbulence to the receiving hydrophone accounted 
for a O.lW signal loss and scintillations due to surface 
swell accounted for 0.07$ signal modulation. The effects 
of bubbles were small because the measurements were rela- 
tively deep and the acoustic frequency was much lower than 
the peak bubble resonant frequency. Measured loss over the 
2m path was on the order of 5$ and phase shifts of up to 
30° occurred for these worst-case conditions. 

Current sound propagation models based upon empirical 
data have undoubtedly included these effects in surface 
scattering or reverberation. But the roughness description 
for the sea surface used by sonar models is based upon sea 
state. The losses here described are, instead, related to 
the swell wave-induced water particle velocities and internal 
waves . 
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APPENDIX A 



Spectral Analysis Plots 

This appendix includes examples of cross-spectral 
analysis plots so that qualitative comparisons between 
sound, temperature and waves during the experiment can be 
made. Included are spectra for 1) sound amplitude versus 
temperature, 2) sound phase versus temperature, 3) wave- 
related effects and *0 comparison of one sound channel 
with another for the same run. The plots included are 
listed below. 



Sound Amplitude Versus Temperature 



Run # Channels Compared 

1 SAMI - TMH2 
1 SAM2 - TMH2 
1 SAM3 - TMH2 



2 SAMI - TMH2 (see text) 

2 SAM2 - TMH2 

2 SAM3 - TMH2 



12 


SAMI - 


TMH2 


12 


SAM2 - 


TMH2 


12 


SAM3 - 


TMH2 



Sound Phase Versus Temperature 

Run # Channels Compared 

1 SPH1 - TMH2 



2 

2 



SPH1 - TMH2 (see text) 
SPH2 - TMH2 



12 



SPH1 - TMH2 
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Wave-Related Effects 



Run # 

1 

1 

2 

6 

12 



Channels 

TMH2 - 
SPH1 - 

SAMI - 

TMH2 - 

TMH2 - 



Sound Channel Comparisons 



Run i? 



Channels 



2 

2 

2 

2 



SAMI - 
SAMI - 
SAM2 - 
SPH1 - 



Compared 



WAVE 

WAVE 

WAVE 

FVBY 

WAVE 



Compared 



SAM2 

SAM3 

SAM3 

SPH2 
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APPENDIX B 



Analyses of Variance Regressions 

Variances of 56 sec samples of data which included 1*100 
data points each were calculated along 26.67 min segments 
of certain channels of runs one, two and 12. Regression 
analyses for sound amplitude and phase versus temperature 
were computed using these variances. The resultant plots 
of the variances with the regression curves are shown after 
Table B5. A summary of the regression equations is given 
in Table Bl. The following mean variance equations were 
computed from the average of the regression equations listed 
for runs one and two: 



V 2 (SAM1) = 3.2 x 10 
V 2 (SAM2 ) = 1.6 x 10 
V 2 (SAM3) = *1.1 x 10 



-3 


v 2 (t) 


(Bl) 


-3 


v 2 (t) 


(B2 ) 


-*J 


v 2 (t) 


(B3) 



p 

If V (T) is eliminated from the above equations (stationary 
assumption), then because SAMI and SAM2 are 30 cm apart and 
60 cm separates SAM2 and SAM3: 



V 2 (SAM) cC (|) L/3 ° cm 
where L = sound path length. 



(b*d 
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This equation can be checked by comparing the sound 
amplitude variances for specific 56 sec samples (where the 
sample is stationary) and similarly deriving the relation- 
ship 



V 2 (SAM) c* (^-r L/x (B5) 

2 

where x is the distance required to reduce V (SAM) by 0.5. 
This is done in Tables B3 - B5 . The following equations 
result : 



? V-L/34.5 

Run 1: v (SAM ) cc ± cm (B6) 

Run 2: V 2 (SAM)<^ |" L/i|2 cm (B?) 

Conversion from base ^ to base e and placing the path 
length reference at lm vice 2m yields (from equations B1 , 
B4, and B5): 

V 2 (SAM) = (0. 3^pb 2 /°c 2 )V 2 (T)e _L/x (B8) 



where : 

L = acoustic path length in meters, 
x = temperature spatial scale for which 
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V 2 (SAM) = V 2 (SAM) x=o /e 
2 

V (T) = temperature variance 

p 

V (SAM) = variance of sound amplitude. 

And the spatial scales in (BD , (B6) and (B7) become *l3cm, 

50 cm and 6lcm, respectively. 

2 2 
Since V (SAM) is known to be a function of a , the 

_3 

variance of the index of refraction, and a = 2.1 x 10 
for run two (which had the best regression fits): 

V 2 (SAM) = 7.7 x 10 ^ a 2 V 2 (T) e _L/x (B9) 



Assuming that x << L, the following approximation can be 
made : 



e 




(BIO) 



Substituting Equation B(10) into Equation B(9) and multi- 
2 

plying by L to compensate for the spherical spreading loss 
yields : 



V 2 (SAM) = 7.7 x 10 2 * a 2 V 2 (T) xL 



(Bll) 
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which is similar to the exponential form of the Born 
approximation: 



v 2 , 2 2 T 

V = k a aL 
o 



(B12 ) 



if k Q is similarly included in Equation (Bll): 

V 2 (SAM) = 1.1 x 10“ 5 k Q 2 a 2 V 2 (T) xL. (B13) 



The plots following Table B5 are included to allow 
qualitative comparisons of conditions between runs and of 
the variances between channels during the same run. For 
each run analyzed the variance versus time plots are given 
followed by regression plots between various channels. 
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TABLE Bl. Regression Equations for Sound Versus Temperature 
Variances 



AMPLITUDE 



Channel 



Run# 



Equation 



SAMI 


1 


V 2 (SAM1) 


= 


C3.^V 2 (T)] X 10“ 3 




2 


v 2 (sami) 


= 


[3.0V 2 (T)] X 10~ 3 




12 


v 2 (sami) 


= 


[l] _ 2V 2 (T) ] x 10~ H 


SAM2 


1 


V 2 (SAM2 ) 


= 


[2 . OV 2 (T) ] x 10~ 3 




2 


V 2 (SAM2 ) 


= 


[1. 2V 2 (T) ] x 10" 3 




12 


V 2 ( SAM2 ) 


= 


[0.7 + 1. 5V 2 (T) ] x 10 


SAM3 


1 


V 2 (SAM3) 


= 


[M . 8V 2 (T) ] x lO -4 




2 


V 2 (SAM3) 


= 


[3 . 3V 2 (T) ] x lO"* 




12 


V 2 (SAM3) 


= 


[0.2 + *IV 2 (T)] x 10 _!< 



-5 



PHASE 



Channel 

SPH1 



Run# 




Equation 


1 


V 2 (SPH1) 


= 1 + illV 2 (T) 


2 


v 2 (sphi) 


= 0.2 + 0. 7V 2 (T) 


12 


v 2 (sphi) 


= 3.7 + 77V 2 (T) 



SPH2 



1 

2 

12 



V 2 (SPH2) 
V 2 (SPH2) 
V 2 (SPH2 ) 



2 + 8V 2 (T) * 
2 + 22V 2 (T) 
2 - 25V 2 (T) 



*From 20 minute record. 



103 



TABLE B2. Run 1 Variances 



sample 


Run 1 


Run 1-18 


Run 1-13 


Run 1-25 


variance 


26 2/3 min 


56 sec 


56 sec 


56 sec 


v t 2 (°c 2 ) 


2.2 x 10" 2 


16.1 x 10 -2 


3.2 x 10~ 2 


0 . *4 3 x 10" 2 


V SAMl^ b ^ 


11.0 x 10 -5 


52.0 x 10“5 


12.0 x 10“5 


18 x 10“5 


v s«<" b2) 


7.6 x 10 -5 


39.0 x 10“5 


7.2 x 10"5 


13 x 10“5 


V SM3^ b2 ^ 


2.5 x 10-5 


9 x IQ" 5 


3.3 x IQ -5 


3.2 x 10"5 


TABLE B3. 


Sound Amplitude 


Distance Loss Factors 


(Run 1) 


sanple 


Run 1* 


Run 1-1 


Run 1-2 


Run 1-3 


distance 


26 2/3 min 


56 sec 


56 sec 


56 sec 


30 cm 


1. *45 


1.33 


1.67 


1.38 


60 cm 


3.0*4 


4.33 


2.18 


*4.06 


90 cm 


*4. *40 


5.78 


3.6*4 


5.63 



quantity 

distance 


mean 

loss 

factors 


standard 
deviation 
in loss 
factors ' 


normalized 

loss 

factors 


accepted 

loss 

factors 


30 cm 


l.*46 


0.18 


1.00 


1.15 


60 cm 


3.52 


1.17 


2. *41 


2 . 30 


90 cm 


5.02 


1.19 


3 . -4*4 


3. 


*not included 


in mean 


loss factors 







10*4 



TABLE B4. Run 2 Variances 



Sanple 


Run 2 


Run 2-9 


Run 2-26 


Run 2-20 


variance 


26 2/3 min 


56 sec 


56 sec 


56 sec 


CM 

O 

o 


46.2 x 10" 2 


98.1 x 10" 2 


143.4 x 10" 2 


121.4 x 10" 2 


V SAMl^ b2 ^ 


140.x 10 -5 


311 x 10" 5 


422 x 10" 5 


420 x 10" 5 


V SAM2^ yb ^ 


56.0 x 10 -5 


122 x 10" 5 


146 x 10" 5 


178 x 10" 5 


V SAM3^ yt)2 ^ 


17.0 x 10" 5 


30.5 x 10 -5 


40.3 x 10~ 5 


46D x 10“ 5 


TABLE B5. 


Sound Amplitude 


Loss Factors with Distance (Run 2 


sanple 


Run 2* 


Run 2-9 


Run 2-26 


Run 2-20 


distance 


2 6 2/3 min 


56 sec 


56 sec 


56 sec 


30 cm 


2.5 


2.5 


2.9 


2.4 


60 cm 


3.29 


4.0 


3.62 


3.9 


90 cm 


8.24 


10.2 


10.47 


9.1 



quantity 

distance 


mean 

loss 

factors 


standard 
deviation 
in loss' 
factors 


normalized 

loss 

factors 


accepted 

loss 

factors 


30 cm 


2.6 


0.3 


1.0 


1.4 


60 cm 


3.8 


0.20 


1.5 


2.8 


90 cm 


9.9 


0.73 


3.8 


4.2 


*not included 


in mean 


loss factors 


• 
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UHRIRNCE FIT 
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032 



033 



235 



027 



032 



023 



Run) ix 



X-SCRLE..--1.Q0F-O1 UNITS INCH* = V 2 
V-SCRLC--5.00C -06 UNITS INCH* - V* 

REGRESSION PLOT 



(TMHX) (V) ^ 

' { SrtMl ^ ( M6LT^> ) 
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'J£G C £0 S 2 C G 2 Q STC GTO 




RUM IX 



X-SC.RLf.--l.00r -01 UNITS INCH- 
y-SCRL.f--S.00f -0? UNITS INCH. 

REGRESSION PLOT 



\i z ( rMHi) (°c z ) 

V Z (S/}M2) (^VOLTS') 
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Rum ix 



K-SRRU DOr-O.l UNITS INCH- * V 2 ( (°c x ) 

y-SCRir--5. OOF -06 UNITS INCH. = v l (s^3) (mol-tv 1 ) 

REGRESSION PLOT 
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RUis) IX 



X-SCPilX'-l. 00T-01 UNITS INCH. = v 1 Ctmhi') (°t*) 
V-SCRLO-l • 00T.+01 UNITS INCH--- vMSPHl) U«v««& 

REGRESSION PLOT 



) 
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K-SCRLf'-l.OOf-Ol UN] TS INCH. » v'Itmhi) (V) 
V-5CALf>S. 00E-02 UNITS INCH- = vMsPHa) (de 5 r«e; i ) 

RtGRhSS I ON PLOT 



1^0 






00 ? 



0 , 



Run) 1 2- 



ft-SCRLF -•! - OOf-Ol UNITS INCH* = ^ 
v-sr,RL F*--*5. npf-02 UNITS INCH.= ^ 

REGRESSION PLOT' 



( w fiVE.) ( 
(TMH2) ( e c"') 
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APPENDIX C 



Spatial Correlation and Structure Functions 

The values of the spatial correlation function and 
structure function calculated during spatial stochastic 
analyses are tabulated in Tables Cl - C6. The functions 
were calculated along the vertical thermistor array. 

The variance of the temperature at each thermistor over 
the ensemble is also given. 

The Tables include: 



Table £ 
Cl 
C2 
C3 
C4 
C5 
C6 



Ensemble 

Run 1, 26.67 min 

Run 2, 26.67 min 

Run 3> 26.67 min 

Low turbulence intensity 

Middle turbulence intensity 

High turbulence intensity 



TABLE Cl. Run 1, Vertical Spatial Correlation and 
Structure Function. 



Thermistor No. 

1 

2 

3 

4 

5 

6 
7 


Variance 

0.5999 
1-0970 
0.0540 
0. 3090 
0.2722 
0.1359 

O.0747 


(°r) 


Correlation 

Function 


Structure 

Function 


Lag 


1.0000 


0.0 


0.0 


0.9635 


0.2081 


5.0000 


0.3745 


0.5841 


15.0000 


0.3654 


0.4673 


20.0000 


0.5507 


0.3978 


25.0000 


2.6243 


O.986O 


25.0000 


0.5208 


0.5926 


30.0000 


0. 1736 


0.3503 


45.0000 


0.2716 


0.2302 


55.0000 


0.2668 


0.7407 


80.0000 


0.2430 


0.9232 


85.0000 


0.0721 


0.3287 


100.0000 


0.5209 


3.1869 


135.0000 


0.9274 


3.3103 


150.0000 


0. 9724 


3.0852 


155.0000 


0.3622 


1.6778 


160.0000 


0.6725 


1.9936 


175.0000 


0.6874 


1.8289 


180.0000 


0.5094 


3.5094 


180.0000 


0.3640 


1.9750 


205.0000 


0.1944 


3.9140 


235.0000 


0.1514 


2.1753 


260.0000 
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TABLE C2. Run 2, Vertical Spatial Correlation and 
Structure Function 



Thermistor No. 



Variance (°C 2 ) 



1 

2 

3 

4 

5 

6 

7 



0.2154 

0.1750 

0.3544 

0.5783 

0.3910 

0.4338 

0.5170 



Correlation Structure Lag 

Function Function 



1.0000 


0.0 


0.0000 


1.0592 


0.3607 


5.0000 


1.0090 


0.3676 


15.0000 


0.9452 


0.0160 


20.0000 


1.0354 


0.0388 


25.0000 


1.0891 


0.4103 


30.0000 


0.9706 


0.0749 


45.0000 


1.1503 


0.1313 


55.0000 


1.0532 


0.2159 


80.0000 


1.0844 


0.6326 


85.0000 


0.9813 


0.2662 


100.0000 


0.5463 


0.2614 


135.0000 


0.6164 


0.6939 


150.0000 


0.5559 


0.3357 


155.0000 


0.5010 


0.4553 


160.0000 


0.5551 


0.9197 


175.0000 


0.4985 


0.5539 


180.0000 


0.5580 


0.4413 


180.0000 


0.4878 


0.6848 


205.0000 


0.5458 


O .6783 


235.0000 


0.4557 


0.9619 


260.0000 



TABLE C 3 . Run 12, Vertical Spatial Correlation and 
Structure Function 



Thermistor No. 



Variance (°C^ ) 



1 

2 

3 

4 

5 

6 

7 



0.2792 

0.2433 

0.1690 

0.1444 

0.1275 

0.0743 

0.0083 



Correlation Structure Lag 

Function Function 



1.0000 


0.0 


0.0 


1.0558 


0.0211 


5.0000 


1.1642 


0.1301 


15.0000 


1.0859 


0.1540 


20.0000 


0.7096 


0.1640 


25.0000 


1.5349 


1.0283 


25.0000 


0.7463 


0.2231 


30.0000 


0.7552 


0.3984 


45.0000 


0.1447 


0.3588 


55.0000 


0.1852 


0.6946 


80.0000 


0.1969 


0.8039 


85.0000 


0.1995 


0.9920 


100.0000 


0. 8057 


1.6228 


135.0000 


0.6778 


1.6854 


150.0000 


0.6295 


1.6494 


155.0000 


0.6143 


2.2873 


160.0000 


0.5148 


2.2931 


175.0000 


0.4773 


2.2355 


180.0000 


0.4121 


I.6676 


180.0000 


0.2998 


2.1737 


205.0000 


0.1487 


1.6762 


235.0000 


0.1214 


2.0123 


260.0000 



TABLE C4 . Low-Intensity, Vertical Spatial Correlation 
and Structure Function 



2 

Thermistor No. Variance (°C ) 



1 

2 

3 

4 

5 

6 
7 

Correlation 

Function 


0 . 0666 
0.0774 
0.0062 
O.OO 38 
0.0028 
0.0114 
0.0014 

Structure 

Function 


Lag 


1.0000 


0.0 


0.0 


1.0785 


0.2173 


5.0000 


1.3019 


0.9860 


15.0000 


1.0129 


1.1896 


20.0000 


0.2793 


4.5337 


25.0000 


9.4721 


32.6599 


25.0000 


0.2274 


5.0118 


30.0000 


0.0791 


6.1497 


45.0000 


0.1556 


4.2863 


55.0000 


0.4158 


0.6737 


80.0000 


0.4850 


0.9095 


85.0000 


0.6269 


1.4669 


100.0000 


3. 1748 


23.5975 


135.0000 


1.6307 


25.8445 


150.0000 


1.5337 


25.6640 


155.0000 


2.9846 


20.1207 


160.0000 


1.5724 


22.1038 


175.0000 


1.1432 


22.5880 


180.0000 


0.7332 


30.3575 


180.0000 


-1.6798 


31.3264 


205.0000 


0.5584 


27.1199 


235.0000 


1.1856 


22.0083 


’ 260.0000 
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TABLE C5. Mid-Intensity, Vertical Spatial Correlation 
and Structure Function 



Thermistor No. 

' 1 

2 

3 

4 

5 

6 
7 



Variance (°C^) 

0.3603 x 10“^ 
0.3726 x 10”f 
0.8536 x 10 1 
0.1329 
0.1011 

0.8996 x 10 | 
0. 3400 x 10" 1 



Correlation 

Function 


Structure 

Function 


Lag 


1.0000 


0.0 


0.0 


1.03^5 


0.2462 


5.0000 


0.9265 


0.3068 


15.0000 


0.8702 


0.1042 


20.0000 


0.8345 


0.2212 


25.0000 


0.0002 


0.0034 


25.0000 


0.8742 


0.4570 


30.0000 


0.6962 


0.3425 


45.0000 


0.2549 


0.7168 


55.0000 


0.2131 


0.9103 


80.0000 


0.2374 


1.1769 


85.0000 


0.1709 


0.8393 


100.0000 


0.0130 


0.8224 


135.0000 


0.0205 


•1.2779 


150.0000 


0.0134 


0.9769 


155.0000 


-0.0008 


0.8464 


160.0000 


-0.0005 


1.3162 


175.0000 


- 0.0008 


1.0017 


180.0000 


0.0149 


0.8641 


180.0000 


-0.0004 


0.8911 


205.0000 


0.0154 


0.3093 


235.0000 


0.0005 


0.3356 


260.0000 
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TABLE C6. High-Intensity, Vertical Spatial Correlation 
and Structure Function 



O 

Thermistor No, Variance (°C ) 



1 

2 

3 

5 

6 
7 

Correlation 

Function 


0.1800 x 10 
0.2462 x 10“ 
0.4578 x 10“ 
0.1231 x 10 
0.1226 x 10“ 
0.5298 x 10“ 
0.2655 x 10" 

Structure 

Function 


5 

3 

3 

1 

2 

2 

1 

Lag 


1.0000 


0.0 


0.0 


0.6971 


9.6513 


5.0000 


0.3813 


9.6565 


15.0000 


0.5125 


0. 3485 


20.0000 


1.7985 


1.7245 


25.0000 


0.0147 


0.1730 


25.0000 


1.1782 


12.0107 


30.0000 


0.8521 


2.9908 


45.0000 


8.3755 


9.2314 


55.0000 


3.7301 


15.2006 


80.0000 


2.9532 


25.7984 


85.0000 


1.9051 


18.2240 


100.0000 


0.2184 


0.1375 


135.0000 


0.2522 


9. 7420 


150.0000 


0.2551 


0.6906 


155.0000 


0.0156 


0. 3438 


160.0000 


0.0156 


10.0158 


175.0000 


0.0170 


0.9675 


180.0000 


0. 3564 


3.8096 


180.0000 


0.0179 


4.2871 


205.0000 


1.0541 


19.7533 


235.0000 


0.0543 


21.5538 


260.0000 
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